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As seen in recent experiments the temperature 
power spectrum for thermally driven turbulence shows a 
strong dependence on the Prandtl number Pr = vj k. v 
is the viscosity of the fluid and k the thermal conductivity 
(diffusivity) of the advected scalar. While for the helium 
celi (Pr « 0.7) there was a large scaling range ||], no 
universal scaling could be found for water [Pr « 7) 0| . 

We took these experiments as a motivation to examine 
the Pr dependence of self similarity features of a passive 
scalar field 0(x, t). The passive scalar d could be the tem- 
perature or a dye (then Pr is often denoted as Schmidt 
number S'è), convected by a turbulent, isotropie veloci ty 
field u(x, t), 



(1) 



/e (x, i) is a forcing term replacing the boundary condi- 
tions. 

Self similarity in turbulence is commonly characterized 
by the scaling exponents of power spectra or of structure 
functions. Alternatively, one may consider the fractal 
dimensions ò^f^ of the d-dimensional graphs of hydrody- 
namic fields. For the passive scalar 0(x, t) the scaling 
exponents of the structure functions are defined by 

(|0(x + r,i)-0(x,i)r') cxr^™\ (2) 

The 

scale dependence of the Hausdorff volume 'tì^'^\G(B'f^y) 
of the graph G{Bf') = |(x, 6i)|x G B^\Q = 0(x)| over 

a ball B^f^ of radius r defines the fractal dimension ò^g^ , 



B'^''\G{Bf^))a^r'o 

A3Ì 



(3) 



In particular, Sg is the Hausdorff dimension of the 



passive scalar graph over a three dimensionai ball Br 



(3) 



x(l) 



The fractal dimension 5g^^ of turbulent passive scalar sig- 
nals is calculated from the fluid dynamical equation. 5g^^ 
depends on the scale. For small Prandtl (or Schmidt) num- 
ber Pr < 10~^ one gets two ranges, 5g^^ — 1 for small scale 
r and 5g^'=5/3 for large r, both as expected. But for large 
Pr > 1 one gets a third, intermediate range in which the sig- nected by 
nai is extremely wrinkled and has = 2. In that range 
the passive scalar structure function Dg{r) has a plateau. We 
calculate the Pr-dependence of the crossovers. Comparison 
with a numerical reduced wave vector set calculation gives 
good agreement with our predictions. 



?(3) 



is the Hausdorff dimension of 



a turbulent signal in space (for fixed time) or - by the 

Taylor hypothesis M - in time for fixed position. 

(tì) 

For the passive scalar the scaling exponent Ci of the 
structure function and the fractal dimension Si'^^ are con- 



<d+(i-Cf^ 



(4) 



We suppose as in |^ that the inequality is in fact sharp. 

In this paper we shall calculate the fractal dimension 
Sg^^ and thus via (||) also the scaling exponent ([^^ from 
the dynamical equation (|^). The main tool of our calcu- 
lation is the volume formula for the passive scalar graphs 
Q , which was introduced as a very useful tool info fluid 
dynamics by Constantin and Procaccia Its main 

advantage is that one can apply rigorous techniques with 
controlied approximations, based on the fluid dynami- 
cal equations. By extending the calculations of |3[g] we 
are able to go beyond an estimate of the exponents: we 
handle also the amplitudes. Doing so, we confirm the 
various scaling ranges addressed in Q and calculate the 
Pr dependence of the crossovers. Comparison with ex- 
periments, simulations, and former theories is discussed. 

For convenience, we measure the passive scalar field 
in multiples of its rms, 9 — 9/9rms- We do not need to 
assume that upper bounds 9max or Umax for the passive 
scalar or velocity field exist, loosening thus the assump- 
tions made in We only need >C2-integrability which 
anyhow is necessary for the existance of structure func- 
tions. 

According to geometrie measure theory [^,|| the Haus- 
dorff volume i7(3)(G(pf ^)) of the graph rè(x) over the 



ball B; 



(3) 



Br reads 



H^-'\GiBr)) 



72fl|2 



(5) 



In one dimension eq.(^) is the well known formula for 
the length of the curve r9{x). Dividing by F^'^-'(Pr) = 
Anr^ /3 and applying Cauchy-Schwarz's inequality we get 



H'^''\G{Br))IV^^\Br)< ^1 + ^ j d%|V2^|2. (6) 

Now the nice idea of |^ was to calculate \V9\^ from 
the heat transfer equation (|l|). In the stationary case it 
is 



|V^|2 = Ì-{-u-V + «;V2}i 



he 



K9r 



(7) 
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We insert the three terms of in (|^) and denote the re- 
sulting terms under the square root in (H) by li, I2, and 
I3 , respectively. Directly from the definition of the ther- 

mal intensity dissipation rate eg, we get It, — — . 

Here 77 = v^^"^ jt^l^ is the Kolmogorov length and e the 
energy dissipation rate /a can further be estimated 
by ^3 ^ PrKe^^l'^ir lr\)^ . Similarily, li can be bound by 

For sufficiently large Re, I2 and can be neglected 
for ali relevant r, i.e., for r smaller than the outer length 
scale L. We are left with li. Applying Gauss' theorem, 
we get 



SiTKr 



^2(x)u(x) • n(x)dA(x). 



(8a) 



n(x) is the unit vector normal to the sphere, directed 
inwards. The centrai idea of now was to introduce 
velocity and scalar field differences in (|8^) to connect 
the Hausdorff dimension of the passive scalar with the 
r-scaling exponents of velocity and scalar differences. 

Here this is achieved in a slightly different way by 
adding a terni ex u(xo) to the rhs of (|8a|), where Xq is 
the center of Br- On average (u(xo)) = 0. Thus we are 
allowed to write 



dB,. 



92(x)(u(x)-u(xo))-n(x)^, (8b) 



where Ar = Anr^ is the surface of the sphere. We again 
apply Cauchy-Schwarz and get 

h < (((u(x) -u(xo)) •n(x))^ 

(8c) 

(. . ■)dBr denotes the averaging over the sphere. The first 
factor under the square root is the flatness of the scalar 
field, which is known to be 3 from experiment ||]. The 
second factor is the longitudinal velocity structure func- 
tion D\\ (r). 

D\\{r) is related to the velocity structure function 
D(r) = (|u(x + r) — u(x)p) via incompressibility 0], 



D\\{r) 



p^D{p)dp. 



(9) 



If we measure lengths in multiples of = u^/'^ / e^/'^ , and 
velocities in multiples of the Kolmogorov velocity w,, = 
[yeY/'^ 1^, D = D/v'^, r = r/rj, we finally obtain 



Vi^KBr) ^ ''''''''' -V 



3V3 , 



Prf^l D\\{f), 
(10) 



which is our main result. The inequalities arise from 
the Cauchy-Schwarz estimation. It is thus reasonable to 



assume that at least the r- and Pr-scaling behaviour is 
correctly given by (M). We thus bave a controlied bound 
for the volume of the passive scalar graphs, 

including the amplitudes. 

Now from experiment it is known that the Batchclor 
interpolation formula is an excellent fit for the velocity 
structure hmction [Bl, 



Dir) 



3(H-a2f2)2/3' 



11.2. 



(11) 



We use the interpolation formula (^l|) to calculate the 
fractal dimension S^^"^ from (fol) and (0) and get 



,(3) 




Prf\ D\\ (f). 



(12) 



The numcrical result from (|T2|) for ólp = — 2 is 
given for several Pr numbers in Fig.l. In principle, four 
cases are possible which we now want to discuss. 

For "small" Pr there are two ranges. If r is sufficiently 
small, we have S'^g^ = 1 and from (|) Cj;^-' = 1. On these 
small scales both the scalar and the velocity field are 
smooth, Z?|| (f) is given by D\\ [f) ~ f^/lS. For increasing 
r a crossover occurs in D\\ (r) and for r > lArj the longitu- 
dinal structure function is given by D\\{r) — (3/ll)fef'»2 ' 

with § ò = a-'^/^l?, = 8.4, Ca"^ near 2/3. The velocity 
field is now fractal but, as Pr is considerd to be small, no 
change can be observcd in the fractal dimension 5^g^ — 1 
of the passive scalar field which stays to be smooth, as 
the 1 under the square root in (10) is stili dominant. 
Physically this means that the diffusivity k is so large 
that the passive scalar field is smooth even on turbulent 
scales. For sufficiently large scale, 




— Pr-^l^ = 0.36Pr-3/4 ^^3^ 



the second term in (|l0| ) becomes dominant. Then (5, 



3/2+C^"V4 and from (^Cr = ^l'^-Ci'l^- Without the 
at most tiny intermittency corrections we have b'^g'' =5/3 
and cj^-* = 1/3. Here both velocity and passive scalar 
field are fractal and scale alike. The classical Obukhov- 
Corrsin scaling theory |l^ is recovered. 

To observe the transition ( |l3|) we must have 
0.36Pr~'^/^ > 14, which implies the condition Pr < 
10^2. Thus "small" Pr are those with Pr < Pri « 10"^. 
Furthermore, if Pr is even smaller than (4^6)"""^, the 
classical Obukhov-Corrsin scaling range can never be 
achieved, because r < L, i.e., the temperature signal is 
smooth on ali scales. 

For "large" Pr the situation is alike for sufRciently 
small and sufficiently large r. But for intermediate r the 
second term under the square root in ( p^ can become 
already dominant although stili r < 14?7, i.e., stili in the 



(1) 
g 
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viscous range of the velocity field, if only Pr is large 
enough. In this range Jg^"* = 2 and consequently cj^^ = 0. 
This means that the passive scalar signal is highly wrin- 
kled although the velocity field is completely smooth on 
that scales. In that situation the dye or the heat is very 
efficiently mixed by the velocity field which is advected 
by the larger turbulent eddies. But since the difFusion 
is very slow (large Pr means small k), concentration or 
temperature differences cannot be smeared out. Simi- 
lar phenomena are obtained when non-turbulent, viscous 
fluids are mixed: fractal pattcrns dcvelop . 

From (|l^) we calculate that this intermediate range 
begins at 



r/v < \/2V5/3Pr-^/^ = 1.22Pr" 



1/2 



(14) 



It ends at r/rj — 14. Thus, to develop such an interme- 
diate range of say, a decade, we must have Pr > Pr^ ~ 1. 
Such Pr we denote as "large". 

Let US now have a look on the passive scalar structure 
fimction Dg{r) = {\9{x + r) — 0(x)p) ex r^2 Neglecting 
again possible intermittency corrections, we have (,2 



2^[^^ = 4 — 'iSy . Fot small r and large r we have C2'' 
and Q =2/3, respectively, as for the velocity structure 
function. But, as derived above, for intermediate r in 
the case of large Pr we have a plateau in the passive 
scalar structure function, Q — 0, Dg{r) = const. The 
complete structure function can easily be reconstructed 
from the scale dependent scaling exponent ('2^'' {''')— 2 ~ 
Sg^\r). The result is shown is Fig.2. 

The plateau has already been predicted by a mean 
field thcory |l2j. In that theory it is observed also 
for Pr > Pru with a very similar Pr^, see Fig.5 of 
The extension of the intermediate range here, 
1.22Pr^^/^ < r/77 < 14, is somewhat difFerent from 
hAQPr^^/^ < r/rj < 15.6Pr^/^ injlll, probably due to 
the mean field approximations in |1^ 

The theory of Batchelor |T3|j9[| also predicts an inter- 
mediate range for Pr > 1. In that theory Dg{r) de- 
pends only logarithmically on r in the intermediate range 
ri < r < r2, but with ri oc Pr~^^^ and r2 independent 
of Pr as in our theory. 

For small Pr < Pri ali three theories predict the tran- 



(1) 



sition from ([^'^ = 1 to Ci"'' = 1/3 for a scale oc Pr 
see eq.(|l|), refs. [||, and ||. 

Note that a plateau in the temperature structure func- 
tion means a fc~^-behavior in the temperature spectrum. 
This fc^^-behavior was also postulated by Kraichnan 

Experiments and full numerical simulations for large 
Reynolds numbers Re and large Prandtl numbers Pr 
are very rare. There are some hints that there in fact 
is a plateau, see the collection of experimental data in 
To get independent confirmation of our predictions, 
we numerically solved eq. (|l|) together with the Navier- 
Stokes equation in a reduced wave vector set calcula- 
tion for several Pr. For details concerning the method, 



-3/4 



see refs. [^5|,Q. The result is shown in Fig.3. For 
small wavevector p, i.e., large r, there is classical scal- 



mg, Q 



C2"'' = 2/3, for ali scales. But for large 



Pr > 1 a plateau develops for medium scales as predicted 
by our theory. For small scales (large p) the spectra fall 
off exponentially, which reflects that the signal is smooth 
{5g^^ = 1) on these scales. 
It becomes particularly evident from our figures, that 



and (■. 



2^-' — A — 2(5g^-' should not be considered as 
global scaling exponent, but, instead, as locai scaling ex- 
ponents Sg^"^ (r) and ('2^'' (r), a concept which we have also 
introduced to examine the intermittency corrections in 
the Navier-Stokes dynamics |16|. If the r-ranges, where 

(9) . 

Q (r) stays at a certain value, are small, it will be dif- 
ficult to extract scaling exponents from experimental or 
simulated data. 

Here this situation appears for Pr > Pru ~ 1. 
This might explain one aspect of the above mentioned 
Rayleigh-Benard experiments of Libchaber and cowork- 
ers where no global scaling exponent can be found for 
the "large" Pr number Pr w 7 [|l|. Note that our the- 
ory can also be applied for an active scalar as for exam- 
ple the temperature field in Rayleigh-Benard convection. 
Eq.(|lO|) remains valid for an active scalar, but, instead 
of the Batchelor interpolation ([ili), the velocity struc- 
ture function of a thermally driven velocity field has to 
be used, which unfortunately is stili not known reliably. 
If we assume Bolgiano-Obukhov scaling ||l^ for the ve- 
locity field. Dir) oc r^/^, we also get Bolgiano-Obukhov 
scaling for the temperature field, Dgir) oc r^/^, if only r 
is sufficiently large. 

The existance of non-global scaling ranges is also pre- 
dicted for low Prandtl number Rayleigh-Benard convec- 
tion, but only for very high Rayleigh numbers |l8| . 

After having introduced global scaling as a paradigma 
in nonlinear dynamics, we might have to get used to non- 
global scaling ranges as the normal case, where univer- 
sal exponents can hardly be derived from experiments 
or simulations, as the individuai scaling ranges are too 
small. 
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FIG. 1. Scale dependent fractal dimension ^g^-* of the pas- 



sive scalar signal for several Pr. 
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FIG. 2. Temperature structure func- 

tion De{f) = De{f)l{eet~^''^v'^l'^). The plateau can easily 
be recognized for sufRciently large Pr. log De (r) is obtained 
by numerical integration of i^j"' {f) over log f . 

FIG. 3. Temperatur spectra from an approximate solution 
of the dynamical equations (by a reduced wave vector set 
method) for several Pr. The Reynolds number always is 
Re = 2 ■ 10*. The dashed-dotted line denotes the velocity 
spectrum of the advecting fluid. fc~'>-bchavior in the calcu- 
lation here with discrete wave vectors means fe~^~^-behavior 
in a continuous wave vector calculation. 
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